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Dispersion  of  Sonnet  in  • Turbulent  Flow 

V.  A.  Krasilnikov  and  V.  I.  Tatarsky 


Dispersion  of  sound  in  a turbulent  flow  is  caused  by  the  fluctuations  of 
the  velocity  of  the  flow  and  by  the  fluctuations  of  the  temperature,  distrib- 
uted irregularly  in  the  flow.  We  shall  assume  that  a turbulent  flow  can  be 
described  by  the  equations  of  an  incompressible  viscous  fluid.  In  an  incom- 
pressible fluid  the  dispersion  by  temperature  irregularities  and  the  dis- 
persion by  fluctuations  of  the  velocity  of  the  flow  are  independent,  The 
total  effect  of  dispersion  is  the  sum  of  these  effects, 

Sound  propagation  in  a turbulent  flow  is  described  by  the  wave  equation 
of  the  acoustics  of  a moving  medium: 


where  f is  the  potential  of  the  sound  field,  c the  Speed  of  sound,  and  v the 
velocity  of  the  notion  of  the  mediian.  In  a system  of  coordinates  that  moves 
with  the  sVerage  velocity  of  the  flow,  v refers  to  the  turbulent  velocity. 

If  we  assune  that  the  fluctuations  of  the  velocity  of  the  flow  are  consider- 
ably smaller  than  the  speed  of  sound  in  the  medium,  we  can  disregard  the  term 
containing  v in  Eq.  (1).  This  equation  then  becomes 
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In  this  expression  and  hereafter  we  omit  the  sumaation  sign  from  1 to  3 
for  the  indices  that  are  repeated  twice.  In  order  to  compare  the  values  of 
the  first  and  second  terms  in  the  right-hand  side  of  (2),  we  must  exmrine 
the  composition  of  the  spectrum  of  the  velocity  field  of  the  flow.  If  the  high- 
est frequency  contained  in  the  spectrum  of  v is  Q,  then  dvi/dt  is  of  the 
order  Qv and  dV/dt  is  of  the  order  of  UP,  where"  u is  the  sound  frequency. 
Thus,  if  the  condition  u » Q is  satisfied,  we  can  disregard  the  second  term 
in  the  right-hand  side  of  Eq.  (2).  According  to  experimental  data,  the  value 
of  Q does  not  exceed  100  cycles/sec,  and  the  second  term  in  the  right-hand 
side  of  Eq.  (2)  will  have  an  effect  only  for  the  lowest  sound  frequencies. 

(It  must  be  noted  that  for  the  frequencies  comparable  with  the  frequencies 
of  turbulent  fluctuations  the  use  of  Eq.  (1)  is  justified  only  sben  the 
fluctuations  of  pressure  in  the  flew  considerably  exceed  the  soad  pressures.) 

The  equation 

JT-Jjr  =7rt» 

(without  the  term  containing  the  acceleration)  can  be  solved  by  the  method 
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of  successive  approximations.  As  a first  approximation,  one  chooses  the  plane 
monochromatic  mare 

with  the  amplitude  A0,  the  frequency  <•>,  and  the  wave  vector  h,-.  Obukhov1 
(assuming  the  statistical  homogeneity  and  iaotropiam  of  the  flow)  obtained 
the  following  equation  for  the  mean  square  of  the  anplitude  of  the  dispersion 
in  the  direction  of  the  single  wave  Vector  a: 


(2SW*)1  l 
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where  V is  the  magnitude  of  the  dispersed  volume,  r the  distance  from  this 
volume  to  the  reference  point,  and 

/?iy  fr)  = «i(r)t>y(r  + p).  (4) 

We  define  the  z axis  by  the  vector  k - fcn  and  the  x,  z plane  by  the  vectors 
k and  k - kn.  If  we  denote  by  0 the  angle  between  k and  a (the  angle  between 
the  directions  of  the  incident  and  dispersed  waves ) , we  obtain 

k — An  «=  |o,  0,  2Asin-|-j;  k = jAcos-|-.  0,  *sin-|-}  ; 

p = p (co8<p  sln<J>,  sin  sin  cos  <|»}  = pm.  ^ 

Here  p,  <p,  and  y are  the  spherical  coordinates  of  the  point  p.  In  terms  of 
these  coordinates,  Eq.  (3)  becomes 

Bit  as  we  know  (aee,  for  instance,  reference  2), 

Rij(p)  =*  (Rrr  — Ru)  tmtnj  — Rttlij,  (7) 

where jRrr(p)  and  Rtt(p)  are  functions  only  of  the  distnce  between  the  points 
in  question.  In  order  to  evaluate  (6)  we  need  only  the  components  flU(  fli8, 

R„ a,  since  only  kt  and  ka  are  different  from  zero.  They  have  the  form 

Ru  (p)  — (Rrr  — Rtf)  cos*  ? sin*<|»  + Ru;  /?!,  (p)  = (Rrr  — Ru ) cos  ?8in<|i  cos 

RU  (P)  - (Rrr  - Rtt)  COS*  <j>  + Ru.  (8) 

In  an  incaapressible  fluid  Rrr  and  Rtt  are  related  by  Karman’ a equation, 
namely 

Ru  “ Rrr  + Vi fRrf  (9) 


If  we  use  this  relation,  we  can  establish  the  form  of  the  dependence  of  jfi  | 8 
on  6 (the  indicatrix  of  dispersion)  without  defining  the  form  of  the  function 
Hrr(p).  When  we  introduce  the  notation 
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If  the  functions  Ar-(p)  an 

p throughout  the  volume  V,  then  the  integration  over  p in 
free)  0 to  « (thie  requires  that  the  condition  V > I*  fee  satisfied;  here  1 is 
the  distance  at  which  the  correlation  function  decreases  considerably). 

•Since  die  contains  the  factor  cos  9,  it  is  abvioss  that  lie  * 0.  As  we 
Shall  nos  show,  Eq.  (9)  implies  that  !*■  ie  alao  sero.  Ssbstituiag  in  (10) 
the  dss  from  (8)  and  integrating  over  9 and  f,  we  find  that 
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where  a * 2fc  sin  0/2.  Wien  we  replace  d«t(p)  by  drr'(p)  + 
grate  by  parts  the  tense  containing  R'rr(p),  we  find  that  Jas 
p^rr(p)*a 

Thus,  Its  * Jos  * 0,  and  hence  ill)  be  cases 


(p)  and  inte  - 
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We  Shall  show  now  that  fu(0)  * Itl(it)  * 0.  for  this  purpose,  let  us  inte-. 
grate  the  expression  for  Jn  (0)  ” Itl  (n)  ever  9 and  W and  again  use  Eq.  (9). 

We  find  that  f|i(0)  B Ih'Cr:)  * 0 if  only  llmpV?rr(p)  —=0. 

- SOD 

Thus,  it  follosM  Iran  the  incompressibility  Fes  Cures  Of  the  twtouleat 
notion  of  a flaid  that  the  indicatrix  of  sound  dispersion  has  srinime  for 

0 0 and  0 ■ It 

for  further  oaaputations,  we . Wat  detemine  the  ions  of  the  faction  Rrr  %>). 
We  assase  for  it  'the  expression 

(14) 

Here  «*  is  the  neen-equare  value  of  the  turbalcnce  velocity  of  the  flow,  and 

1 the  scale  of  the  correlation  defining  the  average  nagnitade  of  the  floc- 
tnetions.  We  can  find  R.t*(o)  fron  Eo.  (9).  Snbstithtiag  thane  eapreowioaa  in 
(10),  'we  find  Jn(0).  For  | | * we  obtain  (oantting  the  inteoeediate  calcs  - 
lations) 


As  exanplos,  Figs.  1 and  2 shew  two  iadioatrioea  of  disparates  correspond 
ing  to  the  valoes  *<1  » 1/^2  mA  M ■ 10.  4haa  hi  »l,  alanst  all  dispersion 
is  directed  forward;  the  -nswiaan  dispersion  is  ia  the  dtaactian  % * l#5kl 
snd  ’the  tendency  toward  narrower  range  of  diaper sian  iscwesaa  nary  rapidly 
with  the  fraqaeacy.  Die  effective  range  of  dispersion  ia  the  aalid  angle  dQ 
'per  wait  of  the  distance  traversed  'by  the  incidmt  ’wave  ia  Jo  uni  Iran  fib.  (15): 


Integrating  (16)  over  all  the  values  of  the  angles,  we  find  the  dispersion 
coefficient  2d: 


Hence,  for  kl  »)  we  obtain 


(18) 


Since  at  high  frequencies  the  main  part  of  the  diaperaed  energy  is  directed 
forward,  we  may  be  interested  here  in  the  part  directed  backward.  Let  us  find 
the  corresponding  dispersion  coefficient,  which  we  designate  2a.  When  kl  >1, 
we  have 
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Thus  the  coefficient  of  the  "backward  dispersion’1  becomes  constant  at 
high  frequencies.  When  we  are  interested  in  the  phase  relations  in  a Wave, 


Fig.  1.  hi  « \ffi.  Fig.  2.  kl  • 10. 

we  aust  take  the  value  of  2d  as  the  dispersion  coefficient  (since  phase  re- 
lations vary  in  a dispersed  wave,  depending  upon  the  magnitude  of  the  wave 
vector,  etc. ).  However,  if  the  given  wave  interests  us  only  from  the  view- 
point of  energy,  we  must  take  as  the  dispersion  coefficient  a considerably 
manlier  value,  namely  2a. 

We  note  here  that  if  the  tens  .allowing  for  the  acceleration  of  the  flow 
is  kept  in  the  original  Eq.  (2),  then  the  following  extra  term  appears  in 
the  expression  (IS)  for  the  mean  square  amplitude  of  the  dispersed  wave: 

V f w»  \ (1  — cos  8)»  (20) 

where  w*  is  the  mew-square  value  of  the  fluctuations  of  acceleration  and 
lt  is  the  correlation  scale  of  the  field  of  accelerations.  In  deriving  (20), 
the  correlation  function  of  the  field  of  accelerations  was  taken  in  the 

form  fly(p)  =w/(r)w/(r+p)=  (a„—  a«)  immj  + a»iu  end  aM  (p)  = V* 
while  arr  was  detenined*  from  the  equation  arr  ■ att  + palt. 


s 


It  is  quite  possible  that  the  dispersion  of  sound  by  the  field  of  accel- 
erStions  may  serve  to  explain  better  s dasping  of  inf resound  in  the  atmos- 
phere, which  damping  is  greater  than  that  accounted  for  by  the  theory  of  dis- 
persion through  viscosity  and  thermal  conductivity. 
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